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We consider the time evolution of a two level system (a two level atom or a qubit) and show 
that it is characterized by a local (in time) gauge invariant evolution equation. The covariant 
derivative operator is constructed and related to the free energy. We show that the gauge invariant 
characterization of the time evolution of the two level system is analogous to the birefringence 
phenomenon in optics. The relation with Berry-like and Anandan-Aharonov phase is pointed out. 
Finally, we discuss entropy, environment effects and the distance in projective Hilbert space between 
two level states in their evolution. 

I. INTRODUCTION 

The role played by gauge fields in quantum theories is of crucial physical relevance, to the point that the gauge 
theory conceptual and formal scheme has become the gauge theory paradigm in quantum field theory (QFT) of 
condensed matter physics and high energy physics. Quantum fluctuations characterize the dynamics of quantum 
systems. Due to quantum fluctuations, quantum systems live on a variety of microscopic configurations. However, 
a remarkable stability is often observed at the mesoscopic and macroscopic scale. The well known QFT strategy to 
face the problem of the emergence of mesoscopic and macroscopic stability out of the fluctuating configurations of the 
system at the microscopic level consists in requiring the local gauge invariance of the system lagrangian. In practice, 
one requires that terms proportional to 9 M 0(x, t), arising from the kinetic term in the Lagrangian when the elementary 
component field 0(x, i) undergoes the local phase transformation </>(x, t) — > </>'(x, t) — cxp(ig#(x, t))</)(x, i), must be 
compensated by the transformation of the gauge field A M (x, t) — > A'^(x, t) — (^(x, t). In such a sense, the gauge field 
thus behaves as a 'reservoir' for the <ft field system [J-Q. This is a well known story, indeed. 

The problem of the stability of macroscopic complex systems arising from fluctuating quantum components is of 
special interest when one considers, for example, the process of defect formation in the non-equilibrium symmetry 
breaking phase transitions characterized by an order parameter Q (e.g. vortices in superconductors and superfluids, 
magnetic domain walls in fcrromagncts, monopoles and cosmic strings in high energy physics and cosmology (5l4ll|. 
the onset of phase locking among the electromagnetic modes and the matter components in the formation of coherent 
domains Q, etc.). In these cases, the gauge structure of the theory reveals itself to be essential and the topological 
characterization of the macroscopically behaving defect appears to be a global feature arising from the quantum 
fluctuating system components. 

In this paper we show that the gauge theory paradigm, in conjunction with geometric phase properties, may be 
applied also to the time evolution of a two level system such as a two level atom or a qubit, thus being relevant also 
in quantum optics and quantum computation studies. As a matter of fact, it has already been realized fl2| that the 
geometric (Berry-like) phase and the related (gauge field) connection pH ] play a relevant role in quantum computing. 
Although such a role of the gauge field has already been recognized, our aim in this paper is to provide, through the 
explicit construction of the covariant derivative, for simplicity, in the case of a two level system or qubit, the physical 
implications of such a gauge structure by showing how it is related to thermodynamical operators such as the free 
energy operator, and how it provides the analogy of the two level system with the birefringence phenomenon in optics. 
The novel picture of time evolution of the two level system thus emerges from our discussion: the two level system 
appears to be embedded in a gauge field background and it evolves in time in such a way that invariance under local 
in time gauge transformations is preserved and it simulates the propagation through a birefringent medium. Finally, 
we also compute the static and dynamic entropy and the distance between the two level system or qubit states in the 
Hilbert space in terms of the Fubini Study metric. 
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II. TIME EVOLUTION AND THE GAUGE STRUCTURE 



In our analysis we consider for simplicity the familiar example of a two level system, e.g. a two level atom or a qubit 
or any other system which might be described by the standard orthonormal basis of two unit (pure state) vectors |0) 
and |1), — Sy, i,j = 0, 1. They may be thought as being cigenstates of the operator 

H =wi|0>(0|+w a |l)(l| , (1) 

with eigenvalues wi and w 2 , respectively: H\0) = wi|0) and H\l) = W2|l). We assume that |0) and |1) arc non 
degenerate eigenstates of H, i.e. that uj\ 7^ In other words, u)\ and o>2 denote the (two different) values of the 
quantum number (energy, or charge, or spin, etc.) characterizing the states |0) and |1), respectively. In the following 
we will think of them as the energy (the frequencies, in natural units h = 1 = c) eigenvalues. 

By a convenient rotation in the plane {|0), |1)}, one may then prepare, at some initial time io = 0, the superposition 
of states 



|0) = a\0) + (2) 
H>) = -p\0) + a\l). (3) 

As usual, orthonormality requires that the coefficients a and f3 satisfy the relations \a\ 2 + |/3| 2 = 1 and a*(5 — af3* = 0. 
Thus we may set in full generality a = e ni cosO and /3 = e' 72 s'm9, with 71 = 72 + mr, n = 0, 1, 2.... Notice that \(f>) 
and are not eigenstates of H due to the fact that u>i ^ u>2- 

In general, in the preparation process we have a limited control (or even no control) on the fluctuations of the a 
and {3 coefficients (the initialization problem [l4Tfl6l ]). However, in some cases, such as in nuclear magnetic resonance 
and electron spin resonance systems, a good precision may be reached in the control of the initialization problem (l7| . 

At time t we have: 

\4>{t)) = e-' lHt \(f>(0)) =e- Wlt (cos9\0) +e-^ 2 - ull}t sme\l)) , (4) 
\^(t)} = e- ifft |V(0)) =e- i " ,lt (-sin0|O) + e-^-^cosflll)), (5) 

where we have used the notation 1 0(0)) = \4>) and | -0(0)) = and, for simplicity, we have considered real a and /3 
(71 = = 72). Of course, time evolution preserves the orthonormalization of the states \4>) and \ip) at any time t. 
By inverting Eqs. (jH) and we obtain the expression for H at any t: 

h =w^\<j,{t)){ ( j>{t)\ +f>wW)){i>(t)\ +w^d0(t)xv(t)i + \m)m)\) , (6) 

where w^^, uj^ and ui^ are given by the time- independent expectation values 

= wi cos 2 6» + w 2 sin 2 6» = (<f>(t)\ id t \<j>(t)) , (7) 

ui^ = u)! sin 2 6» + w 2 cos 2 6> = (ip(t)\ id t \tp(t)) , (8) 

uo H , = i(wa - wi) sin 29 = (i/>(t)\ id t \</>(t)) , (9) 

and ui^ = ui^fp. Note that the oj^ "mixed term", responsible for "oscillations" between the states \4>(t)} and \ip(t)), 
appears in H since w 2 — &i 7^ and the "mixing angle" 9 is non- vanishing 1 . 

We thus have H\(p{t)) = uj^ \4>{t)) + ui^ \tp(t)) and H\ij}(t)) = uj^, \i(>(t)) + ui^ 1 0(i)) . On the other hand, since 
H\cj>{t)) =id t \<t>{t)) and H\if>(i)) =id t \ijj(t)) (cf. Eqs. g]), © and (|T])), we get the evolution equations 

id t \((t)) = u d \C(t)) + !£(*)>> (10) 

where \((t)} denotes the vector doublet |C(*)) = (I^W) , {^(t)}) 1 " , tod = diag(Lo^ : uj^) and <7\ denotes the Pauli 
matrix. By using the notation g = tan 29 = U J^' , with 6u> = ui^,^ — w<^, we have u/^ = ^gSui. We also put 
Aq = Aq 1 ' (7 1 = i 5uj o\. Then we may write 

D t = d t + iuj H ,ai = d t + igA^ai, (II) 



1 Same situation occurs in the mixing of neutrinos and in general of particles with different masses [T3l . 
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which acts as the covariant derivative, where g and play the role of the coupling constant and the (non-abelian) 
gauge field, respectively (a similar situation occurs in the different context of neutrino mixing, see ref. [l9j). The 
motion equations (fTU|) now can be written as 

iA|C(*)> = w d |C(t)>. (12) 

iD' t \(>(t)) = w d \C'(t)), (13) 



It is easy to show that 
with 



D' t = d t + ig{A^<n +d t X(t)a 1 ), (14) 

|C'(t))=e-^M*)- 1 |C(i)), (15) 

so that, defining U(t) = e r %9X{t)ai , it is 

U(t)(iD t \((t))) = iD' t U(t)\((t)) (16) 

and 

gA^' ax = U{t)gJ$<TiU-\t) + i{dtU{t))U~\t) , (17) 

as it should be indeed for a gauge field transformation (see Eq. (fT4")l ). 

We can express the above result by saying that the time evolution of the vector doublet \C,{t)) (our two level system 
or qubit) is controlled by its coupling with a non-abelian gauge field background so to preserve the invariance of the 
dynamics against local in time gauge transformations (phase fluctuations). 

We also note that since the only non- vanishing component of A^ is A$ and this is a constant (Aq = ^Suiax), 
the field strength F^ v is identically zero. This is a feature which, for example, occurs in the case where the gauge 
potential is a pure gauge (with non-singular gauge functions). 



A. The gauge field background as a birefringence medium 

We now show that the time evolution described above can be interpreted in terms of a birefringence phenomenon 2 . 
Contrarily to what we have assumed above, let us assume now that the states |0) and |1) are degenerate states, 
namely their time evolution "in the vacuum" is given by 



|0(t)> \ = / e"*" \ / [0(0)) \ 
\l(t)) ) \ e"-* ) \ |1(0)) ) 1 



(18) 



where u> = 2 tt i>, and the propagation speed in the vacuum is vo = X v. Suppose then that the propagation occurs in 
a medium presenting different refraction indexes, n\ and n 2 for |0) and |1), respectively, i.e. where the propagation 
over a given path of length I occurs in different times, t\ and ti for |0) and |1), respectively: 

i £m £ £n 2 

ti = — = = tni ; t 2 = — = = tn 2 , (19) 

vi v v 2 v 

where v\ and v 2 are the propagation speeds in the medium for |0) and |1), respectively, and t = Time evolution 

is then described by the phase factors e~ tutl = e~' luJlt and e _lw * 2 = e~' luJ2t for the two states, respectively, where 
uiti = uj—rii ~ 2tt v trii = 2-k v.it = u>it, i = 1,2, has been used, together with Xi v = Vi, XiVi = and 
rii = 321 = Thus, for the mixed states \4>) and \ifj) given by Eqs. ^ and §3§ we have 

/ cos(9 e-^-^sine \ ( 10) \ 
= e- lUlt , (20) 

\ - sin 6 e'^-^ cose /VI 1 )/ 

which is the time evolution generated by H given by Eq. ([1]) with cji ^ oj 2 (cf. Eqs. ^ and (|5|) ). In conclusion, Eq. 
([20]) shows that, provided that u)\ ^ oj 2 , for 9 ^ j + the effect of time evolution through the refractive medium is 

equivalent to the effect of the background gauge field A^ = \{u 2 — Wi) cos 29 = }^uj(n 2 — ti\) cos 29, which indeed 
disappears when propagation occurs in the vacuum, ri\ = n 2 = uq = 1 (i.e. u>i = lu = u> 2 ). 




2 The analogy with birefringence has been considered for the case of neutrino mixing |2(I 
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B. The free energy, the topological phases and the Fubini— Study metric 



Let us now write Eqs. (flQ|) as 

(H - w^trOKC*)) = Wd|C(*)>, (21) 

with the covariant derivative denoted by H — uj^, o\ . Then the operator 

F = (H - w^tri) (22) 

may be interpreted as the free energy operator, provided that one identifies the term cu^ o\ = gAo with the entropy 
term TS in the traditional free energy expression, where the "temperature" is T = g and the entropy S = Aq. Wc 
thus see that time evolution is controlled by the free energy (f2"2")l where the gauge field plays the role of the entropy. 
In terms of the states \4>{t)) and \ip{t)), the term TS is written as: 

ts = Lo H ,{\4>{t))(m\ + mwm ■ (23) 

In order to better understand this feature, it is convenient to consider the geometric phases associated to the system 
evolution. We observe that it is an easy matter to compute the Berry-like phase. Indeed, one immediately gets the 
geometrical phase (3$ [2l] |: 

P+ = <P+ f ((f>(t)\id t \<f>(t))dt = 2nsm 2 9, (24) 



which is independent of the u>i's, i = 1,2, and depends only on the "mixing angle" 9. In Eq. ([24| . p = — 2 ™^ x , 
u>i 7^ ll>2, and we exploited the fact that after a period T = LJ 2 * Ui it is \4>(T)) = e lv \(f>(0)). 
Similarly, for the state \i^(t)) one finds (3^, = 27rcos 2 9 and thus (3^ + fitf, = 2ir for any 9. 

Another geometric invariant is the Anandan-Aharonov phase discussed in ref. [22[ . It has the advantage to be well 
defined also for systems with non-cyclic evolution and is given by s = 2 J Auj(t)dt, where Aui(t) is the variance given 
by 

auj 2 = a^ = Aw^ = m\H 2 \m) - mwm)) 2 = ^ 2 H , = ^ . (25) 

with £ = <f>, ip . The relation between the entropy and the geometric invariant s is obtained by considering that 

(C(t)|TS<n|C(i)) dt = J (0)\gA^\C(t)) dt = 2 J u H ,dt = s. (26) 

It is interesting to note that the relation between TS and the variance of the energy Alo = ui^ is through the 
non-diagonal elements of H, namely it is proportional to the energy gap, w^ — ^i, between the two levels (cf. Eq. (|25j) 
and ©)• We also recognize that the integrand (((t)\g A^\((t)} in Eq. (|2"6"]l is related to the adiabatic connection R3I 
emerging in the study of the non-abelian holonomy (generalized Berry phase) fl2l ]. Moreover, one can also show |l8| 
that these connections are related with the parallel transport of the vectors in the parameter space, as well known 
[l~2T ] . In the Appendix A we compute the static and the dynamic entropy associated to the reduced density matrix. 

Finally, it is instructive to analyze the invariant s in Eq. (|26j) in terms of the distance between states in the Hilbert 
space. Let us gcncrically denote by either \<fi(t)) or \ip(t)). Their evolution is governed by the Schrodinger 

equation 

i ft|f(t)) = H\£(t)) , t = M. (27) 
Expanding the state \£(t + dt)) up to the second order in dt and taking into account that Jj-ff = 0, we have 

dt 2 

mm + ^)) = i-idt mmm) — (m&im + o(dt 3 ) , 

(28) 

and |(f(£)|C(* + dt))\ 2 = 1 - d^Auj 2 ^ + 0(dt 3 ) . That is 

uj 2 

\{£{tM(t + dt))\ 2 = 1- dt 2 — sin 2 29 + 0(dt 3 ), £ = <t>,1>, (29) 
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where we have used Eqs.([7]), (|5J) and (|25|) and 

m)\H 

(ip(t)\H 2 \tp(t)) = u 2 2 cos 2 (9 + uf sin* 8. (31) 



(<f>(t)\ H 2 \(f>(t)) = uj\ cos 2 9 + lo^ sin2 9 , ( 30 ) 

r2 \ a i,fj.\\ , ,2 f\ i , , 2 



We also have 

[<0(*)|^(t + dt)>l a = |<V(*)|^(* + dt)>I a =d* 2 Aw5 V) + 0(dt 3 ). (32) 

The Fubini-Study metric [22| is defined as 

ds 2 = 2g^dZ»dZ v = 4(1 - |(£(£)|£(i + di))| 2 ), (33) 

where Z M are coordinates in the projective Hilbert space V, which is the set of rays of the Hilbert space T~L. From 
Eqs. (|29| . (|32[) and (|33|). we have the infinitesimal geodetic distance between the points U(\<fi(t)}) and II(|0(£ + dt))) 
in the space V 

ds = 2Au) Si dt = u_ sin26dt. (34) 

The rate of change of this distance is 4f = w_ sin 29 — 2cj0^,, with w_ = UJ2 — lo\ ^ 0. In the case of the above 
two level or qubit states, the Fubini-Study metric coincides with the usual metric on a sphere of unitary radius: 
ds 2 = dO 2 + sin 2 Qdp 2 , with = 2 (6 = mixing angle) and O £ [0,7r]. Since 9 is constant, we have ds = sin 20 dip 
and, by comparison with Eq. p4j) , dtp = w_ dt . We thus obtain 



sin 2$ e?iy2 = 2 J ui^dt, (35) 

which is the Anandan-Aharonov invariant (cf. Eq. (|26|) ). Thus, the Anandan-Aharonov invariant s represents the 
distance between evolution states, as measured by the Fubini-Study metric, in the projective Hilbert space V. 



III. CONCLUSIONS 



In conclusion we have shown that time evolution of a two level system or qubit is controlled by a covariant derivative 
accounting for the coupling of the state with a (non-abelian) gauge field background so to preserve the invariance of 
the dynamics against local in time gauge transformations. We have shown that the effect of the gauge field background 
can be depicted as the effect of a birefringence phenomenon, the gauge field bacground acting as the analogous of the 
refractive medium. We have also shown that the covariant derivative plays the role of the free energy with the gauge 
field acting as the entropy. In such a picture time evolution is controlled by the free energy. Finally, the relation 
of our result with the geometric phase and the so-called adiabatic connection has been pointed out. We have also 
shown that the distance in the projective Hilbert space between two level or qubit evolution states is measured by 
the Fubini-Study metric in terms of the Anandan-Aharonov geometric invariant. For the reader convenience and for 
completeness, further comments on the entropy and the environment effects are reported in the Appendix A. 



Appendix A: Mixed states, entropy and environment effects 

By exploiting the Schmit decomposition theorem (see e.g. 0, 53), one ma y always "double" the system under 
study; denote it by A. The "doubled" system, denoted by A, is introduced in such a way to work in the composite 
Hilbert space H A ^ = %a ® H-^ with states |^ ^) = J2 n V™™ l°™ "«) e A' S n w " = 1- The density matrix for 
mixed states of the system A, = ^2 n w n \a n )(a n \, is obtained by tracing out the system A: 

P A = w n \a n )(a n \ = ^2 y/WnWm | On) (dm | Tr (\d n )(d m \) = Tr^ {p ABA ) , (Al) 

n nm 

where the relation (a m \a n ) = 6 nm has been used. Then, one can show that the "tilde" system A can be interpreted as 
the thermal bath or reservoir for the original system A [111 [26| and the free energy and the entropy can be defined 3 . 



3 Such a construction is equivalent to the GNS construction in the C* -algebra formalism and requires the quantum field theory framework 
25, 26]. 
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The state j) is recognized to be an entangled state of the tilde and non-tilde modes and the entropy provides a 
measure of the entanglement [lH, HtJ . 

We now compute the static (linear) entropy for the qubit states |0(i)) and \tp(t)} given by Eqs. dU and ([5]), 
respectively. One introduces the states |0) and |1) as 

|0) |0>®|6), (A2) 
|1) -+ |1>®|1>, (A3) 

and uses Eqs. (|A2[) and (|A3[) in Eqs. Q and (JS|). The density matrices for the states in H s where S = {0, 1} and 
S = {0, 1}, are denoted by p^ = !£(*),£(*))(£(£), £(i)l where £ = 0, V and £ = 0, V- The reduced density matrix p^ 
(and similarly for p^) is obtained by tracing out the tilde-system S, and vice- versa. Thus one obtains: 

Pf = Tr § P^ = cos 2 9 |0) (0| + sin 2 9 11) (1| , (A4) 
p| = Tr 5 p - cos 2 9 |6) (6| + sin 2 (A5) 

The static linear entropies Sl associated to the reduced matrices p^ and p^ are then: 

Sl14] = 2(1-Tr s [(4) 2 ])=sin 2 (20), (A6) 
S L [pf] = 2(1-Tr5[(p|) 2 ]) = sin 2 (20). (A7) 

Recall that sin 2 29 = ^jj- (cf. Eq. ([9])), where cj_ = W2 — wi 7^ 0. Let us now compute the dynamic entropy. 

Consider the state |0(i)) in Eq.@ (we can proceed in a similar way for \tp(t))) and express it in terms of the states 
10(0)) and |V(0)): 

|0(t)) - A^(t)|0(O)) + A^(i)|V(0)), (A8) 
where e are the amplitudes: 

A H (t) = (0(0)10(0) = e^" 1 * cos 2 9 + e" 8 " 2 * sin 2 9, (A9) 
= <"0(O)|0(t)> = e~ lUlt sin (9 cos 9 + e^ 2 ' sin (9 cos 9, (A10) 

respectively The tilde-states |0) and |V) are introduced, for any t, as : 

|0) -> |0)® |0), (All) 
|V) -> |V)®I^>- (A12) 

The reduced density matrices are now, for any t, 

pf = Tr^ P4> = |A^(t)| 2 |0) (0| + \A^(t)\ 2 |V) (V|, (A13) 
pf = Tr nP4> = |A^(i)| 2 |0) (0| + |^(t)| 2 |V) (VI, (A14) 
where 1Z = {0, tp} and 72. = {0, V}- The dynamic entropies Sl are 

S L (p£) = 2(l-Tr 7e [(^) 2 ])=4|^(0| 2 |^(i)| 2 =4F^(t)P^(<), (A15) 
S L ( P f) = 2(l-Tr^[(pJ) 2 ])=4|^(0| 2 |A^(i)| 2 =4F^(i)P^W, (A16) 
where P^^{t) and P^^{t) are the probabilities of the transitions 0—^0 and — > V: 

P*_^(f) = sin 2 (20)sin 2 (^^-^ , (A17) 



iW*) = 1 - sin 2 (20) sin 2 ( i^_J^t ) , ( A 18) 



respectively. 
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